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Magnetohydrodynamic  Modelling  of 
Interplanetary  Disturbances  Between 

the  Sun  and  Earth 


I.  INTROIM'CTION 


The  roles  of  solar  llat  .-s  an  1  si  i  ••am  -st  ream  interactions  and  their  effects 
upon  th»  solar-  -a  i n;  1  an  1  interplanetn  rv  magnetic  field  (IMFl  are  now  generally 
appreciated  .vithin  both  'he  seien'ifio  community  and  the  larger  group  of  diverse 
operational  users  of  sophisticated  high -technology  svstems.  The  latter  group  is 
especially  concerned  with  magnet ospheric  and  ionospheric  disturbances  during 
geomagnetic  storms.  Thus,  the  need  for  long-range  forecasting  (on  a  time  scale 
of  many  hours  or  even  several  days)  has  become  clear.  An  MUD  model,  together 
with  several  operational  spacecrafts  that  can  ntovide  continuous  "ground  truth" 
in  the  solar  wind,  is  an  obvious  candidate  far  a  key  role  in  providing  this  need. 

This  paper's  purpose  is  to  demonstrate  the  capability,  unique  at  this  time  to 
our  knowledge,  of  a  self-consistent  MUD  model  to  produce  operational  forecasts 
that  are  specifically  directed  to  the  following  questions: 

(1)  Will  a  solar  disturbance  affect  the  geophysical  environment  (specifically 
ionospheric  and  magnetospherie  disturbances)?  ...  and,  if  so, 

(2)  When  will  a  geomagnetic  storm  start0 

(3)  How  long  will  it  last  ° 

(4)  How  severe  will  if  tie" 


(Keceived  for-  publication  December  2f't!2> 


It  k.  emphasized  that  the  MIID  model  that  will  be  discussed  is  a  development  in  the 
long  tradition  of  simulations  of  initial-boundary  value  problems  that  are  solved 
with  numerical  finite  difference  solutions  of  partial  differential  equations.  The 
model  is  not  a  kinematical  approximation  as  performed,  for  example,  by 
Hakamada  and  Akasofu  for  the  kind  of  heliospheric  problem  addressed  here. 

2 

More  specifically,  we  will  use  the  model  devised  by  Wu,  Dryer,  and  Han.  This 
non-planar  model  (so-called  because  it  includes  all  three  components  of  the  solar 

wind  plasma  and  IMF)  is  the  most  recent  development  that  follows  upon  the  ear- 

3  4  5 

lier  work  of  Wu,  Han,  and  Dryer;  Dryer,  Wu,  and  Han;  D'Uston  et  al;  and 

g 

D'Uston.  The  earlier  work  clearly  demonstrated  the  development  of  fast- 
forward  (as  well  as  reverse)  MHD  shock  waves,  their  asymmetrical  propagation, 
as  well  as  their  multiple  interactions  (as  in  the  well-known,  multiple-flare  situ¬ 
ation  in  August  1972).  It  was,  however,  restricted  to  only  two  components  of  the 
plasma  bulk  velocity  and  IMF  in  the  ecliptic  plane  as  a  consequence  of  its  purely 
2-D  formulation.  As  part  of  our  long-range  program  (that  includes  a  fully  3-D, 
time-dependent  formulation),  a  non-planar  extension  to  the  earlier  2-D  studies 
was  incorporated  (Wu,  Dryer,  and  Han).  This  "quasi -3  D"  or"2-l/2D"  formu¬ 
lation  now  makes  it  possible  to  relate  solar  disturbances  to  geomagnetic  disturb- 

7 

ances  via  the  studies  reviewed,  for  example,  by  Smart,  Garrett,  and  Shea.  We 
have  in  mind,  more  specifically,  the  use  of  solar  wind  parameters  as  proposed 

g 

by;  (1)  Perreault  and  Akasofu  who  advocate  the  use  of  the  energy  coupling  function, 

1.  Hakamada,  K.  ,  and  Akasofu,  S. -I.  (1982)  Simulation  of  three-dimensional 

solar  wind  disturbances  and  resulting  geomagnetic  storms.  Space  Sci. 

Rev.  3^:3 -70. 

2.  Wu,  S.T.,  Dryer,  M. ,  and  Han,  S.  M.  (1 982)  Non-planar  MH D  model  for 

solar  flare-generated  disturbances  in  the  heliospheric  equatorial  plane, 
submitted  to  Solar  Physics. 

3.  Wu,  S.T.,  Han,  S.M.,  and  Dryer,  M.  (1979)  Two-dimensional,  time- 

dependent  MHD  description  of  interplanetary  disturbances;  simulation  of 
high  speed  solar  wind  interactions,  Planet.  Space  Sci.  27:255-204. 

4.  Dryer,  M. ,  Wu,  S.  T.  ,  and  Han,  S.M.  (1980)  Two-dimensiona-1,  time- 

dependent  MHD  simulation  of  the  disturbed  solar  wind  due  to  representative 
flare-generated  and  cor  onal  hole-generated  disturbances,  Geofis.  Intern. 

19:1-15. 

5.  D'Uston,  C.  ,  Dryer,  M.  ,  Han,  S.  M. ,  and  Wu,  S.  T.  (1981)  Spatial  structure 

of  flare-associated  perturbations  in  the  solar  wind  simulated  bv  a  two- 
dimensional  numerical  MHD  model,  J.  Geophys.  Res.  80:525-534. 

0.  D'Uston,  C.  (1981)  Dynamiquo  des  Ondes  dc  Choc  Interplanetaires 

Fngendrees  par  des  Kruptions  Solaires,  PhD  thesis  in  Physical  Sc iences, 
University  Paul  Sabatier-  de  Toulouse,  Toulouse,  France. 

7.  Smart,  D.  F.  ,  Garrett,  H.B.,  and  Shea,  M.A.  (1980)  The  prediction  of  A  K, 

ap,  and  Dst  at  time  lags  between  0  and  30  hours,  in  Solar -Ter- res  trial 
Predictions  Proceedings.  Vol.  II,  R.F.  Donnelly,  Fd. ,  U.  S.  Government 
Printing  Office,  Washington,  D.  C.  ,  pp.  399-413. 

8,  Perreault,  P.  ,  and  Akasofu,  S.  -I.  (1978)  A  study  of  geomagnetic  storms, 

Geophys.  J.  R.  Astr.  Soc.  54:547. 


8 


e;  and  (2)  Burton,  McPherron,  and  Russell^  who  advocate  the  half-wave  rectifier 
concept  for  the  magnetosphere  as  represented  by  the  quantity  VBg  [  where  Bg  is 
the  southward  (negative)  IMF  component  ...  corresponding,  later,  to  our  positive 
component,  Bg] .  Other  workers  have  examined  the  statistical  use  of  these  (and 
other)  quantities  for  the  prediction  of  geophysical  parameters  such  as  the  AE, 

Dst,  AL,  and  Kp  indices  (see,  for  example.  Smart,  Garrett,  and  Shea;  Baker 
et  al;^  Clauer  et  al;^  Holzer  and  Slaving  Kamide;^  and  Kamide  and  Akasofu^). 

There  is  no  useful  purpose  served  if  the  present  paper  were  to  advocate  one 
index  in  relation  to  another.  We  view  our  role  as  one  that  will  provide  all  of  the 
parameters  that  are  used  to  compute  the  cited  indices  once  we  are  provided  with 
knowledgeable  (even  educated  guesses!)  conditions  near  the  sun  on  a  regular 
basis.  1  This  latter  phrase  really  means  a  time  scale  that  is  dependent  on  solar 
activity.  ]  The  former  phrase  refers  to  solar  observations  that  can  delineate  not 
only  solar  flares  (radio  sweep  frequency,  soft  X-ray),  but  also  those  that  can 
infer  the  slower  changes  such  as  coronal  hole  (optical,  magnetograph)  develop¬ 
ment  or  demise.  Only  time-dependent  modelling  has  the  capability  for  fulfilling 
this  requirement.  Also,  strictly  speaking,  a  3-D,  time-dependent  MHD  model 
is  necessary  when  it  is  known  that  large  latitudinal  changes  (compared  to  those 
in  the  heliospheric  equatorial  plane)  take  place.  We  will  consider  this  problem 
in  the  future.  Stated  mathematically,  our  present  quasi-3D  model  will,  then, 
assume  that  partial  derivatives  of  all  dependent  variables  (density,  temperature, 
velocity,  IMF)  with  respect  to  the  heliocolatitude,  0,  are  negligible  in  comparison 
with  those  derivatives  in  the  radial  (r)  and  heliolongitudinal  ((A)  directions. 

Studies  in  MHD  modelling  of  the  interplanetary  environment  are  incomplete 
without  direct  confrontation  with  observations.  It  is  essential  that  the  MHD 
models  be  tested  in  relation  to  the  "ground  truth"  provided  by  solar  wind  plasma 
and  magnetometer  instrumentation  aboard  spacecraft  outside  the  influence  of 


9.  Burton,  R.K. ,  McPherron,  R.  L. ,  and  Russell,  C.  T.  (1975)  An  empirical 
relationship  between  interplanetary  conditions  and  Dst,  J.  Geophys.  Res. 
8£:4204.  “ 

10.  Baker,  D.  N. ,  Hones,  Jr.,  E.W.,  Payne,  J.B.,  and  Feldman,  W.C.  (1981) 

A  high  time  resolution  study  of  interplanetary  parameter  correlations  with 
AE,  Geophys.  Res.  Lett.  8:179-182. 

11.  Clauer,  C.  R. ,  McPherron,  R.L.,  Searls,  C. ,  and  Kivelson,  M.G.  (1981) 

Solar  wind  control  of  auroral  zone  geomagnetic  activity,  Geophys.  Res. 
Lett.  8:915.  - 

12.  Holzer,  R.  E. ,  and  Slavin,  J.  A.  (1982)  An  evaluation  of  three  predictors  of 

geomagnetic  activity,  J.  Geophys.  Res.  87:2558-2562. 

13.  Kamide,  Y.  (1982)  Comment  on  "An  evaluation  of  three  predictors  of 

geomagnetic  activity",  J.  Geophys.  Res,  (submitted). 

14.  Kamide,  Y. ,  and  Akasofu,  S.  -I.  (1982)  Notes  on  the  auroral  electrojet 

indices  (in  press). 
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Earth's  bow  shock.  We  have  performed  such  studies  in  the  past  as  summarized 
in  several  review  articles  (Dryer;  ^  ^  Wu^'  and  a  doctoral  thesis  (D'Uston).^ 
Further  studies  with  the  present  quasi-3D  model,  however,  are  still  needed  before 
it  can  be  used  as  an  operational  tool  by  AWS  forecasters.  Fortunately,  multi¬ 
spacecraft  observations  by  the  Helios,  Voyager,  Pioneer,  Prognoz-8,  ISEE-3, 
and  IMP  spacecraft  are  now  being  provided  by  various  experimental  groups 
(Burla^a  et  al;  Dryer  et  al;  H  Intriligator;  Klein  and  Burlaga;  Richter 

-t  al;  Zastenker  et  al;^  Intriligator  and  Miller;^®  and  Borrini  et  al.  These 
groups  are  delineating  properties  such  as  the  global  extent  of  flare-generated 
shock  waves,  energetic  particle  modulation,  IMF  distortions  into  "magnetic" 
clouds,  and  post -shock  plasma  properties  in  everincreasing  detail.  We  plan  to 
make  comparisons  in  the  future  of  multi-spacecraft  observations  with  the  predic¬ 
tions  provided  by  our  MHD  model. 

Prior  to  such  a  massive  undertaking,  however,  it  is  necessary  that  general 
properties  of  some  canonical  examples  be  provided  by  the  model  so  that  major 
physical  properties  can  be  identified.  Such  is  the  purpose  of  the  present  report. 
Therefore,  we  will  examine  the  interplanetary  disturbances  that  are  produced  by 
a  canonical  set  of  flare-generated  shock  waves. 


2.  ANALYSIS 

The  governing  equations  for  the  time-dependent,  non-planar  MHD  problem 
are  mathematical  expressions  of  the  physical  conservation  laws.  For  convenience 
of  numerical  computation,  they  are  written  in  quasi-conservation  form  in  spheri¬ 
cal  coordinates  and  specialized  for  the  equatorial  plane  (0  =  90°).  The  equatorial 
plane  may  be  either  the  ecliptic  plane,  the  solar  equatorial  plane,  or  even  the 
equatorial  plane  in  a  tilted  heliospheric  coordinate  system.  The  latter  system 

may  also  be  a  candidate  for  the  present  analysis  as  a  consequence  of  the  discovery 

3  1 

of  the  solar  current  sheet  (for  example,  Rhodes  and  Smith  ),  the  empirical  study 

32  1 

of  Zhao  and  Hundhausen,  and  the  kinematic  study  of  Hakamada  and  Akasofu. 

A  schematic  sketch  of  heliospheric  space  is  shown  in  Figure  la  with  the  two 

cross-hatched  planes  indicating  the  meridional  plane  and  the  equatorial  plane  of 

our  analysis  Figure  lb.  We  will  consider  a  180°  segment  of  the  latter  plane. 

The  general  equations  for  the  conservation  of  fundamental  physical  quantities 

(neglecting  dissipation,  except  at  shock  waves!  in  mks  units  are  identical  to  those 
2 

given  by  Wu  et  al  and,  for  completeness,  are  summarized  in  Appendix  A. 


Because  of  the  large  number  of  references  cited  above,  they  will  not  be  listed 
here.  See  References,  page  25. 
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Figure  la.  Illustration  of  the  Spherical  Coordinate  System 
and  Heliospheric  Space.  The  equatorial  plane  is  chosen 
for  analysis 
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Figure  lb.  Schematic  Presentation  of  the  Computational  Domain, 
Flare  Location  at  <j>  -  30°  (Note  That  <b  -  0°  at  j  =  It,  and  Posi¬ 
tions  of  Hypothetical  Spacecraft  Monitors 
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The  ambient  steady-state  of  the  interplanetary  medium  is  obtained  by  setting 
3(  ) /Ot  s  0.  The  result,  within  our-  domain  of  supersonic/super-Alfvenic  flow, 
is  a  1-D  flow  with  a  canonical  Archimedean  magnetic -field  topology  in  the  equa¬ 
torial  plane  (Figure  la).  Initial  conditions,  including  small  meridional  velocity 
and  magnetic  field  components,  at  the  inner  boundary  (taken  at  18  ft  where  Kg 
is  the  solar  radius)  are  varied  until  the  desired  representative  values  of  the 
various  plasma  and  field  parameters  are  obtained  at  1  AT.  The  parameters  at 
1.  10  AT  and  at  18  ft  are  given  in  Table  1.  The  computational  domain  can  be 
adjusted  to  suit  the  particular  problem  of  interest.  The  present  code  can  handle 
problems  in  subsonic  as  well  as  in  the  supersonic  flow  regime;  however,  for 
subsonic  flow  problems,  the  boundary  conditions  should  be  modified  to  satisfy 
the  compatibility  relations  and  non-reflective  boundary  conditions.  The  latter 

requirements  have  been  considered  in  several  sola r-oriented  problems  (see,  for 
33  3  4 

example,  Hu  and  Wu;  and  Wu  et  al.  )  but  not.  as  yet,  in  the  present  code. 


Table  1.  Initial  Steady-State  Solar  Wind  Conditions 


Pa  rn  meter 

18  R 

s 

1.  1  AT 

n 

1.8v  io3  cm'3 

10.0  cm'3 

V 

r 

250  km  /s 

3  60  km/s 

ve 

0.  2  km /s 

4.  4  x  10'3  km/s 

Vd 

4.  0  km/s 

0.  5  km/s 

B 

r 

300  nT 

2.3  nT 

1  nT 

4.  4  x  10'2  nT 

-30.  0  nT 

-2.  0  nT 

T 

1.  1  v  10°  °K 

3.  3  x  104  °K 

With  respect  to  gr  id  size,  we  chose  Ar  =  2Rg  and  Adi  =  3°.  To  cover  180° 
of  heliolongitude  and  radial  distance  just  beyond  1  AT,  111  radial  positions  and 

33.  Hu,  Y.Q.,  and  Wu,  S.  T.  (1982)  A  full-implicit-continuous-Kulerian  (FIC K) 

scheme  for  multidimensional,  transient  magnetohydrodynamie  (MHD) 
flows,  J,  Computational  Phys.  (in  press). 

34.  Wu,  S.  T.  ,  Hu,  Y.Q.,  Nakagawa,  Y.  ,  and  Tandberg-Hansscn,  T.  (1982) 

Induced  mass  and  wave  motions  in  the  lower  solar  atmosphere.  I.  effects 
of  shear  motions  on  flux  tubes,  submitted  to  Ast  rophys.  J  ■ 


til  azimuthal  positions  were  required  (Figure  lb).  'mailer  grid  size  would  result 
in  better  resolution  of  the  flow  field  at  the  expense  of  larger  computer  storage. 
Numerical  results  with  the  present  grid  size  showed  relatively  good  accuracy. 

Karh  of  our  sola  r -fin  re-gene  rat  ed  MUD  shock  perturbations  was  introduced 
at  the  lower  boundary  lasting  a  time,  At,  with  a  ramp  input  and  decay  (to  elimi¬ 
nate  numerical  instability)  at  each  siue  of  the  temporal  duration.  The  Rankine- 
Ilugoniot  MUD  shock  jump  conditions  were  separately  solved  for  an  arbitrarily- 
selected  shock  velocity  taking  the  lower  boundary  steady  state  conditions  (Table 
1)  as  the  pre-shock  conditions.  The  computed  post -shock  conditions  were  then 
taken  at  the  center  (j  11)  of  a  24°  longitudinal  extent  in  our  equatorial  plane. 
Thus,  the  shock's  central  meridian  (CM)  is  located  at  <*>  -  .40°,  measured  count¬ 
erclockwise  from  tin-  horizontal  axis  in  Figure  11).  We  assumed  a  sinusoidal 
decay  of  the  post -shock  conditions  to  either  side  of  this  shock  wave  in  order  to 
simulate  the  existence  of  an  asymmetrical  shock  generated  originally  by  our 
’’canonical"  solar  flare.  We  will  examine  global  (that  is,  lb0°  equatorial  plane) 
snapshots  of  magnetic  field  topology,  density,  and  temperature,  as  well  as  the 
time  series  of  the  various  properties  as  "measurer!"  by  a  hypothetical  observer 
at  two  radial  positions,  0.4  AT  and  1  .A l  ,  along  the  flare's  CM  and  at  two  addi¬ 
tional  hidiolongitudin.il  positions  of  .T!°\V  and  !'0'\v  of  the  CM  (Figure  lb). 

Computations  were  performed  for  a  wide  range  of  shock  velocities,  V  ,  at 
the  center  of  the  are  shown  in  Figure  lb  at  ltil!  100,  1000,  1100.  2000,  2000, 
and  2  100  km  sec  *.  The  temporal  duration,  A',  ot  the  simulated  flare  was  taken 
for  most  of  these  shock  velocities  as  1<!00,  4000,  1400,  and  7200  sec.  These 
latter  figures  are  based,  rather  h.atrist  ically,  on  representative  (characteristic) 
time  durations  of  integrate  !  soft  X-ray  flux  observations  by  earth-orbiting 
spacecraft.  We  will  discuss  only  the  case,  V  .  4000  sec  *  (1400  sec),  in  vhis 

report  in  order  to  locus  attention  on  essential  physical  features  of  the  interplan¬ 
etary  disturbanc  e  between  the  sun  and  Karth. 

.4.  GKNKRAI,  RKSI.l.TS 

This  section  will  consider  some  representative,  essential  results  for  the 
case  mentioned  in  Section  2.  A  detailed  examination  of  the  capability  of  the 
numerical  code  to  produce  physically-relevnnt  additional  quantities  is  given  in 
Appendix  I?  for  the  same  ease. 

For  the  purpose  of  illustration,  we  will  discuss  here  only  the  results  gener¬ 
ated  by  a  "flare"  that  generates  an  interplanetary  shock  wave  with  a  velocity  of 
4000  km  see  ^  at  14  solar  radii.  The  temporal  duration  of  this  large  flare  (as 


indicat  i'll  by  a  hypothetical  sof!  X-ray  ;twi’!  vat  ion,  inti-grab  1  aver  tin-  s >l:o  iiak), 
is  taken  to  be  ,yt  =  f>400  see. 

Figure  2  shows  tile  contours  of  total  pressure  projected  onto  the  Solm  euua  - 
'Hi  it,!  plane  tit  t  0,  20.  2.  40.  2.  and  00.  1  hr  after  initiation  of  tin  flare.  The 
contours  are  labeled  by  a  logarithmic  scale  for  Fog  ^  (  2nk  1  •  with  units 

of  lyties  cm  “.  Ihe  classical  A  echini'  dean  spiral  IMF  is  not  shown  at  '  0  hr. 

It  is,  however,  shown  at  subsequent  times,  superimpose)  upon  the  total  pt  •  ssiirc 
contours,  as  the  flare-generated  shock  wave  expands  toward  Faith's  mbit.  We 
call  attention  to  the  local  minimum  of  the  total  pressure  that  occurs  in  'lie  region 
of  the  most  severe  1M  I  kinking  (note  the  S-shnped  reversals).  |r  is  cb  nr  that 
this  11-field  rotation  is  due  to  the  shock's  passage  through  tie  solar  wind.  (Some 
detailed  physics  of  this  feature  are  discussed  in  Appendix  11.  )  y  et,  within  the 
peripheral  areas  one  detects  total  pressures  that  are  higher  than  the  earlier 
un  iisturbed  values  at  the  same  local  ions. 

l'he  basic  parameter,  solar  wind  velocity,  is  graphics  1 1  ,\  shown  in  Figure  3 
as  an  evolving  2-1.)  surface.  Only  the  quadrant  0  »'  ©  «•  5*0°,  is  shown  so  that  the 
viewer,  in  inertial  space,  can  observe  the  disturbance's  major  features.  The 
undisturbed,  steady -state  solar  wind  (homogeneous  in  heliolongit tide)  is  shown 
in  Figure  2(a).  The  "cliff"  at  the  lower  left  corner  of  each  subset  of  this  figure 
simply  represents  a  discont  imiance  of  the  calculation  beyond  1.  1  A1  and,  there¬ 
fore,  should  be  ignored.  (We  plan  to  extend  the  computational  domain  in  the 
future.)  It  is  seen,  in  Figure  3(b),  that  the  central,  strongest  portion  of  the  fast 
forward  shock  is  about  three-fourths  of  the  distance  to  the  earth's  orbit.  Tile 
fast  reverse  shock's  strongest  portion  is  also  seen  in  this  figure  for  t  20  hr. 
(Vote  that  the  large  velocities,  as  indicated  in  the  figure  title,  required  a  scale 
change  from  that  used  in  the  other  parts  of  tile  figure.)  Attcutuation  of  the  shock 
wave  in  the  heliolongitudinal  direction  is  clearly  seen  for  both  fast  forward  and 
reverse  shocks.  Motion  of  the  fast  forward  shock  out  of  the  domain  (except  for  a 
longitudinal  remnant  at  ~o  ;i0-!'0n)  is  also  seen  in  Figure  3(c).  Die  fast  reverse 
shock  passes  out  of  the  domain  at  the  later  times  (~!i(l  hr)  as  noted  in  Figure  3(e) 
and  3(d).  As  the  flow  gradually  returns  to  the  original  steady -stale  condition,  we 
note  several  radially -orient ed,  relatively -higher  velocity  "  ridges"  that  are  clus¬ 
tered  around  the  radially-extended  IMF.  The  ridges  represent  the  loci  of  several 
slow  mode  waves  that  are  the  nonlinear  counterpart  of  the  well-known  "petals"  in 
the  Friedrich's  phase  velocity  diagram.  It  is  recalled  that  tile  slow  mode  propa¬ 
gation  vectors  are  restricted  in  angular  extent  about  the  local  direction  of  the 
magnetic  field.  Figures  3(e)  -  3(g)  show  the  longitudinal  motion,  away  f'om  the 
original  shock's  CM,  toward  both  lower  and  higher  values  of  o. 

Figure  4  shows,  in  a  time  series  formal,  radial  velocity,  density,  temper  - 
nture,  and  IMF  total  magnitude  at  1.0  AF.  Figures  4(a)  -  4(d)  show  tin  time 
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owes  or  > 


e  2.  Total  Pr  essure  (2nkT  +  B  /8*r)  Contours,  Together  with  Superpositic 
at  t  >  0  hr).  Following  Introduction  of  the  Shock  Wave,  Vs  -  3000  km  sec" 
Numbers  are  logio  (dynes  cm'2) 
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Figure  3.  Three-dimensional,  Polar,  View  of  the  Solar  Wind  Velocity  Sur  face 
Within  the  Quadrant  0  *•  <$  <  90°.  Since  the  eomputat ionnl  domain  is  terminated  at 
1.  1  At',  the  velocity  appears  as  zero  outside  this  distance  in  this  presentation  and 
should,  together  with  the  vertical  "cliff",  be  ignored.  The  outer  boundary  condi¬ 
tions  allow  for  continuous  flow  through  it  ns  indicated  by  the  smooth  passage  of, 
for  example,  the  forward  and  reverse  shocks  through  r  1.  1  AF:  (a)  t  0  hr 
(ambient,  steady-state  solar  wind  as  given  in  Table  1),  (b)  t  20  hr  (complete 
structure  of  the  disturbance  is  container!  within  the  computational  domain.  Note 
scale  change  compared  with  (a).  (AV)niax  (118  km  sec-11,  (c)  t  40  hr  (fast 
forward  shock  has  nearly  completed  its  passage  through  the  1.  1  AF  radial  posi¬ 
tion.  (AV)niax  within  the  computational  domain  is  383  km  sec'll,  (d)  t  -  (>0  hr 
(strongest  por  tion  of  the  fast  reverse  shock  is  approaching  the  1.  1  AF  radial  posi¬ 
tion.  (AV)max  29a  km  sec-M,  (,,)  t  go  hr  (strongest  pot  ion  of  the  fast  re¬ 
verse  shock  has  passed  thr  ough  the  1.  1  AF  radial  position.  (AV)m  19(1  km 
sec-1 1,  (f)  t  100  hr  (fast  reverse  shock  has  completed  its  passage  through  the 
1.1  AF  radial  position.  Multiple,  radially-aligned  "ridges"  are  remnants,  via 
slow  mode  wave  propagation,  of  the  original  disturbances.  (AV)ma„  100  km 
sec  “  1 1 ,  and  (g)  t  100  hr  [continuation  and  gradual  attentuation  of  disturbance 
and  return  to  ambient  conditions.  (AV)max  (10  km  sec'll 
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Figure  3.  Three-dimensional,  Polar,  View  of  the  Solar  Wind  Velocity  Surface 
Within  the  Quadrant  0  --  ©  90°.  Since  the  computational  domain  is  terminated  at 

1.  1  AC,  the  velocity  appears  as  zero  outside  this  distance  in  this  presentation  and 
should,  together  with  the  vertical  "cliff,  be  ignored.  The  outer  boundary  condi¬ 
tions  allow  for  continuous  flow  through  it  as  indicated  by  the  smooth  passage  of, 
for  example,  the  forward  and  reverse  shocks  through  r  1.  1  AU:  (a)  t  =  0  hr 
(ambient,  steady-state  solar  wind  as  given  in  Table  1),  (b)  t  =  20  hr  Icomplete 
structure  of  the  disturbance  is  contained  within  the  computational  domain.  Note 
scale  change  compared  with  (a).  (AV)max  -  018  km  sec'll,  (c)  t  40  hr  (fast 
forward  shock  has  nearly  completed  its  passage  through  the  1.  1  AU  radial  posi¬ 
tion.  (AV)max  "'thin  the  computational  domain  is  383  km  sec'M,  (d)  t  =  80  hr 
(strongest  portion  of  the  fast  reverse  shock  is  approaching  thr  1.  1  AU  radial  posi¬ 
tion.  (AV)mnx  29a  km  sec_M,  (e)  t  ~  80  hr  (strongest  portion  of  the  fast  re¬ 
verse  shock  has  passed  through  the  1.  1  AU  radial  position.  (AV)max  -  190  km 
see"*  I,  (f)  t  -  100  hr  [fast  reverse  shock  has  completed  its  passage  through  the 
1.1  AU  radial  position.  Multiple,  radially-aligned  " ridges"  are  remnants,  via 
slow  mode  wave  propagation,  of  the  original  disturbances.  (AV)mnx  100  km 
sec"M,  and  (g)  t  100  hr  (continuation  and  gradual  attentuation  of  disturbance 

and  return  to  ambient  conditions.  (AV)  -  00  km  sec'M  (U'ontd) 

max 
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TEMPERATURE  IN  K  RADIAL  VELOCITY  IN  KM, SEC 


series  at  the  CM  of  the  flare-generated  shock  pulse.  It  is  seen  that  the  fast 
forward  shock  arrives  at  1.  0  AU  at  about  t  =  24  hr  and  a  more  well-developed 
fast  reverse  shock  arrived  at  about  t  =  60  hr. 


MHD  SIMULATION  AT  r  =  1  0  A  U 
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Figure  4.  lime  Series  of  iindial  Solar  Wind  Velocity,  Density,  Temperature, 
and  dotal  IMJ*  Magnitude  at  1. 0  AF.  CM  of  flare-general ed  shock  (a-d) 


18 


4.  GEOPHYSICALLY  RELEVANT  SOLAR  WIND  PARAMETERS 


As  noted  in  the  Introduction,  there  is  currently  an  evaluation  of  various 
combinations  of  solar  wind  parameters  and  their  statistical  relevance  in  regard 
to  various  geomagnetic  indices.  We  have  already  implicitly  answered  the  first 
of  the  four  fundamental  questions  posed  in  the  Introduction.  We  have  done  so,  at 
least,  for  the  case  considered  thus  far:  a  uniform  steady-state  solar  wind  with 
a  canonical  spiral  IMF  that  is  disturbed  by  a  very  fast  shock  at  181<^,  namely 
the  case  of  V  =  3000  km  sec  \  At  -  5400  sec.  Thus,  if  the  earth  were  located 
at  or  within  90  W  of  the  pulse's  CM,  we  would  reply  in  the  affirmative  to  the 
question:  Will  a  disturbance  (presumed  to  be  known  with  some  confidence  close 
to  the  sun)  reach  the  earth0  We  will  study  the  case  of  larger  angles,  both  west 
and  east  of  the  CM  in  the  future. 

'1’he  second  question  is,  if  a  disturbance  will  reach  the  earth's  magneto¬ 
sphere,  when  will  it  begin0  Arrival  of  the  fast-forward  shock  wave  at  the  mag¬ 
netosphere  will  produce  a  storm  sudden  commencement.  (Whether  a  storm  main 
phase  will  follow  is  discussed  later.  )  We  have  intensively  examined  the  ease 
mentioned  in  the  preceding  paragraph.  As  noted  earlier,  we  also  studied  a  wide 
range  of  initial  shock  velocities  and  flare  duration  times.  The  range  of  initial 

velocities  was  500  -?  V  s  3500  km  sec  *  and  3ti00  «  At  «■  7200  sec.  As  noted, 

s 

details  of  the  disturbed  plasma  and  IMF  for  most  of  these  cases  are  given  in 
Appendix  B.  The  upper  limit  for  transit  time  is  that  given  for  the  ease  of  the 
original  steady-state  solar  wind.  Thus,  the  integrated  transit  time  for  a  test 
particle  moving  at  250  km  sec  *  at  18K^  (Table  1)  is  110  hr  in  all  directions 
(since  CM  in  this  limiting  case  has  no  meaning).  A  test  particle  moving  faster 
than  this  speed  will  obviously  arrive  sooner  at  1  AT.  Thus,  a  small  pulse  that 
produces  only  the  three  well-known  MUD  (fast,  slow,  and  Alfven)  modes  will  not 
produce  a  fast  shock  until  the  disturbance  velocity  exceeds  the  fast  mode.  Our 
lowest  Vg  exceeds  this  bound:  we  have  not  studied  (here)  the  case  of  purely  linear, 
or  even  nonlinear,  MUD  wave  modes  and  their  coupling  (although  our  computer 
code  is  fully  capable  of  handling  such  situations).  Thus,  strictly  speaking,  our 
results  for  shock  transit  time,  shown  in  Figures  5(a)  -  5(c),  apply  only  to  the 
fast-forward  shock  wave.  The  curves'  upper  bound  is  simply  the  110-hr  transit 
time  for  the  undisturbed  steady-state  solar  wind.  The  curves  have  been  extrapo¬ 
lated  to  Vy  4500  km  sec  ^  as  part  of  our  polynomial  fits  to  the  individual  points 
found  from  our  more  restricted  study. 

We  made  no  attempt  to  extrapolate  below  ot  above  the  assumed  set  of  "flare 
duration"  times.  At.  Classically,  a  shock  wave  that  results  from  a  delta  function 
input  of  energy  is  called  a  "blast  wave".  (We  leave  the  interpolation  of  this  case, 
At  -  0,  to  the  reader.)  Rigorously,  then,  all  of  our  shocks  are  really 
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Figure  5.  Shock  Transit  Time  to  1.0  ATJ.  The  initial  shock  velocity  is  as¬ 
sumed  to  be  known  at  18RS;  no  allowance  is  made  for  the  transit  time  from 
1.0-  18HS.  The  "flare  duration"  time  is  assumed  to  be  known,  on  a  heu¬ 
ristic  basis  from  the  temporal  profile  of  the  soft  X-ray  flux  (for  example, 
1-8A)  from  a  given  flare:  (a)  CM,  (b)  33°W  of  the  CM,  and  (c)  90°W  of  CM 


"piston-driven"  shocks  because  At  is  finite.  Since  our  choices  for  At  are  small 

relative  to  the  transit  times  for  most  of  our  Vg  choices,  we  would  be  justified  in 

referring  to  these  shocks  as  "blast  waves"  only  after  At  has  bei'n  completed  and 

tile  elapsed  time  is  large  compared  with  At.  As  a  counter  example,  for  the  case 

of  V  ■»  3300  km  sec'1  and  At  >  7200  sec,  and  at  CM,  the  shock  transit  tine 
s 

will  be  of  order  of  10  hr.  If  At  were  of  order  of  3  to  5  hr,  we  should  refer  to 
such  a  shock  as  "piston-driven"  all  the  way  to  1  AU.  In  any  case,  Figure  5  pro¬ 
vides  a  prediction  (neglecting  an  additional,  but  small,  corr  ection  for  the  transit 
time  from  1-1811,.)  for  the  second,  operationally-motivated,  question:  when  will 
the  disturbance  (read:  fast  forward  shock  for  our  study)  reach  earth0  A  final 
caveat  is  important:  Figure  5  pertains  only  to  the  representative  solar  wind 
velocity  profile  given  in  Table  1,  thus  the  transit  times  for  other  cases  would 
require  revision  upward  or  downward! 

The  third  and  fourth  questions  posed  in  the  Introduction  (how  severe0  how 
long0)  are  much  more  difficult  to  answer  for  the  general  case  of  spatially-in- 
homogeneous  photospheric,  coronal  and  soinr  wind  conditions,  at  best,  and  for 
the  temporal  changes  that  are  superimposed  on  such  a  background,  at  worst.  We 

will  return  to  our  "ansatz"  or  "canonical"  example  of  the  case  V  -  3000  km 

- 1  S 
see  (At  -  3400  Sec)  in  order  to  offer  straightforward  answers  to  these  two 

questions.  Our  answers  are  obtained  self-consistently  from  the  foregoing  results 
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within  the  context  of  the  somewhat  controversial  debate  currently  underway  as 
noted  in  the  Introduction.  We  will,  therefore,  discuss  our  calculations  for  the 
energy  flux  parameter,  e,  and  the  electric  field,  VB^.  Note  that  the  former  is 
evaluated  in  the  frame  of  solar-ecliptic  coordinates  (not  solar-magnetospheric 
coordinates).  Note,  also,  that  our  velocity  is  the  resultant  magnitude  of  the 
velocity  (which  coincides  with  other  current  usage).  Our  "B  "  is  really  Bg,  which 
is  positive  when  directed  southward  (bv  virtue  of  our  use  of  standard  spherical 
coordinates).  Since  B.  ,  as  used  by  magnetospherie  physicists,  is  negative  when 
directed  southward  of  tin'  ecliptic,  the  reader  is  asked  to  change  signs  (mentally) 
when  comparing  these  results  with  those  shown  in  magnetospherie  papers  on 
this  subject. 

Figure  (i  shows  the  time  ser  ies  for  e  and  VBfl  at  1.0  AF  for  the  case:  Vs  = 
3000  km  sec  \  At  ~  3400  sec.  If  we  were  to  adopt  the  recipe  proposed,  for 
example,  by  Ilakamada  and  Akasofu,  ^  wo  would  predict  high  latitude  activity,  as 
measured  by  the  All  index,  to  commence  at  t  r  24  hr  (immediately  following  the 
SC)  when  e  rises  to  the  first  threshold  of  10lfi  er  g  sec  We  assume,  for  this 

exercise,  that  the  earth's  magnetosphere  is  located  at  C'M.  The  second  threshold 

10-1  * 
of  10  '  erg  sec  would  be  breached  at  t  ~  2(>  hr  when  the  magnetospherie  ring 

current  would  develop,  as  measured  by  Dsi,  and  the  main  phase  of  the  geomag¬ 
netic  storm  would  begin  with  associated  auroral  zone  motion  and  ionospheric 
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disturbances  at  lower  latitudes.  Interpreting  the  10  erg  sec  threshold  some¬ 
what  liberally,  we  might  suggest  that  the  Dst  would  fluctuate  somewhat  until  t  ~ 

38  hr,  at  which  time  the  low  to  medium  latitudes  would  recover  to  pre-storm 

1 8 

conditions.  The  high  latitude  activity  would  persist  until  e  falls  below  10  erg 
sec  *  at  t  ~  f> 7  hr.  Overall,  then,  one  would  suggest,  on  the  basis  of  this  exer¬ 
cise,  that  the  activity  (when  earth  is  located  at  CM)  would  primarily  be  located 
at  high  latitudes,  with  a  brief  Dst  perturbation.  As  for  the  fourth  question,  a 
t'me  duration  of  37-24  -  23  hr  would  be  offered  as  a  prediction  for  the  storm's 
duration. 

On  the  other  hand,  if  we  were  to  adopt  the  concept  of  the  half-wave  yectifi e r 
proposed,  for  example,  by  Burton  et  al,  we  would  predict  AK  activity  commenc¬ 
ing  (as  in  the  previous  example)  at  t  ~  24  hr  and  ending  al  t  ~  31  hr  when  VBg  in 
Figure  li  is  positive  (negative,  in  the  parlance  of  the  magnetospherie  physicists). 
The  large  negative  values  in  Figure  r>  would,  in  this  recipe,  not  be  relevant  since 
the  i>?  component  (B„  in  our  parlance)  would  be  northward  during  the  interval 
31  -  t  <  fiO  hr. 

Clearly,  these  two  recipes  discussed  would  also  be  used  for  the  most  prob¬ 
able  cases  when  the  earth's  magnetosphere  is  at  locations  other  than  CM.  We 
hope  it  is  clear  to  the  reader  that  location  relative  to  CM  is  not  necessarily  the 
prime  consideration.  The  strength  of  the  perturbation  (velocity)  as  well  as  the 
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orientation  of  the  component  evolve  in  the  earthward  direction  follow  ins  some 
kind  of  change  at  the  sun.  Hence,  it  is  possible  (as  shown  for  the  33°W  position 
in  the  central  panel  of  l-'igure  t>)  that  e  can  be  greater  at  positions  other  than  CM. 


5.  CONCLUDING  KFMAKkS 

We  have  use.J  a  time-dependent,  nonplanar,  2-1)  Mill)  computer  model  to 
simulate  a  number  of  flare-generated  shock  waves  and  their  concomitant  effects 
on  the  interplanetary  medium.  The  example  of  V  3000  km  sec  \  with  an 
assumed  flare  (or,  equivalently,  energy  input)  time  of  :>400  sec,  was  used  as  a 
detailed  illustrative  case.  A  representative  steady-state,  solar  wind  with  axial 
symmetry,  was  assumed  for  simplicity  as  part  of  this  "ansatz"  problem  with 
three  components  of  the  velocity  vector  and  magnetic  field.  Hence,  all  nonlinear 
wave  mode  coupling  is  included  sol f -consistent  1\  . 

The  follow  ing  major  result  evolved  from  the  presentation  of  the  complete 
sets  of  plasma  and  IMF  parameters.  We  have  established  the  fact  that  the  model 
can  be  satisfactorily  used  to  compute  any  form  of  geophysically -relevant  solar 
wind  parameter  (such  as  e  or  VB  )  for  long-range  (~  days)  predictions  of  geomag¬ 
netic  and  ionospheric  activity.  Future  work  must  include  retrospective  tests 
that  employ  bo'll  solar  and  mult  i -spacer  raft  observations  for  specific  intervals 
of  time.  Perhaps  several  weeks  before,  during,  and  after  one  or  more  solar 
flares  and/or  coronal  hole  evolutionary  observations  should  be  a  minimum  tem¬ 
poral  range  that  should  be  used  in  a  lest  made  for  the  Mill)  model. 

Still  lacking  in  the  present  work  is  the  incorporation  of  additional  physical 

phenomena  such  as  multi-fluids,  dissipation  (thermal  conduction,  viscosity,  and 

electrical  resistivity),  and  more  init  ial  -boundary  value  conditions  that  are  more 

realistically  reflective  of  the  changing  mood  of  solar  conditions.  Perhaps  the  two 

most  important  omissions  in  the  present  work  arc  the  incorporation  oT  hybrid 

:( 1  n 

kinetic /fluid  effects  as  discussed  by  Cuperman  '  and  Wu  and  more  realistic 
considerations  of  heat  flux  and  non-Maxwellian  distribution  functions.  Fast,  but 
not  least,  even  the  most  "primitive”  version  of  the  Mill)  model  as  discussed  here, 
must  be  tested  by  direct  confrontations  with  a  number  of  spacVcraft  observations 
during  specific  (that  is,  not  statistical)  intervals  of  solar  interplanetary  activity. 


.'la.  Cuperman,  S.  (10112)  Solar  wind  theory.  Space  Sci.  Key,  (in  press). 
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p  =  2p  I?T  (A  7) 
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The  first  term  in  Equations  (2)  -  (5)  gives  the  continuity  equation  for  consci - 
vation  of  mass.  The  second,  third,  and  fourth  terms  give  the  momentum  equa¬ 
tions  in  the  radial,  meridional,  and  azimuthal  directions,  respectively.  The  fifth, 
sixth,  and  seventh  terms  give  the  induction  equation  in  the  same  three  directions 
(  again,  9(  ) / 9 a  Ol  with  infinite  electric  conductivity  governing  the  coupled  inter¬ 
action  between  the  plasma  motion  and  the  magnetic  field.  Finally,  the  eighth  term 
provides  the  energy-conservation  equation  for  our  one-fluid  system.  This  dynam¬ 
ic  system  is  assumed  to  be  adiabatic  and  the  symbols  have  their  usual  meaning: 
p  is  the  mass  density;  V  -  (Vr>  Vg,  V^)  denotes  the  mass -flow  velocity  vector; 

T  is  the  temperature;  B  =  (B;,,  Bg,  B^)  denotes  the  magnetic  field;  p  -  2pRT  is 
the  gas  pressure,  accounting  fot  equal  electron  and  proton  contributions  via 


31 


Dalton's  law:  M  .  is  the  solar  mass;  W  the  gravitational  potential;  G  the  gravita 
tionnl  acceleration  at  the  surface  of  the  Sun  (r  I-f  );  ■>  the  specific  heat  ratio 
and  jj  is  the  magnetic  permeability  in  a  vacuum.  The  independent  varic 
hies  are  time  (t)  and  the  radial  (r)  ami  .azimuthal  (0)  coordinates  in  this  non- 
planar  (2-D,  3 -component)  system. 
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Appendix  B 

Additional  Results 

\Ve  have  presented  major  results  in  tin-  main  text  for  the  case:  V  3000 

-1  S 

km  sro  ,  At  5400  sec .  In  order  to  examine  the  more  detailed  physics,  for 

completeness  we  shall  now  outline  some  other  interesting  results  in  this  appendix 
for  this  ease.  The  capability  of  the  MHO  model  to  provide  temporal  and  spatial 
details  of  the  disturbed  interplanetary  plasma  and  IMF  will  also  be  demonstrated. 

Figure  B1  shows  the  [IMF  as  projected  into  the  equatorial  plane  at  t  =  0,  20.2, 
40.2,  and  HO.  1  hr.  Recalling  the  location  of  the  shock's  CM  (equivalent,  essen¬ 
tially,  to  the  flare's  CM),  we  note  the  fast  forward  shock-induced  deformation, 
particularly  along  and  east  of  the  CM.  Of  particular  interest  is  the  reversed 
IMF  polarity  (from  outward  to  inward)  within  the  disturbed  flow  at  t  =  40.2  hr. 

We  will  refer  again  to  this  feature  and  the  physical  mechanism  with  which  this 
kinking  is  associated.  As  the  shock  recedes  into  the  far  reaches  (that  is,  beyond 
1  A  l  in  the  present  case)  of  the  equatorial  plane,  the  IMF  gradually  returns  (I  = 
200  hr)  to  its  pre-flare  Archimedean  topology.  The  numbers  on  the  heavier  lines 
are  constants  associated  with  magnetic  vector  potential;  they  have  no  meaning 
other  than  their  usefulness  in  the  identification  (as  time  progresses)  of  the  foot- 
point  of  a  par  ticular  field  line  at  lttii^,.  These  field  lines  have  already  been 
shown,  in  Figure  2.  superimposed  upon  the  total  pressure  contours.  Here  in 
Figure  Bl,  however,  the  distorted  IMP'  can  be  seen  more  clearly. 

Figure  B2  shows  the  contours  of  I.og^T  together  with  the  superimposed  IMF 
from  the  previous  figure.  The  formation  of  the  fast  forward  MHD  shock  wave  is 
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Figure  Bl.  Interplanetary  Magnetic -field  Line  Deformation  Following  Introduction 
of  a  Solar  Flare-generated  Shook  Wave,  Vs  =  3000  km  sec"*  and  At  =  5400  sec. 

(a)  t  -  0  hr,  (b)  t  20.2  hr,  (c)  t  -  40.2  hr,  and  (d)  t  =  60 .  1  hr.  The  numbers  on 
the  heavier  lines  are  chosen  by  the  computer  (from  the  vector  potential  constanis) 
and  indicate  the  original  longitudinal  position  at  18RS 


now  observed  more  clearly.  We  can  also  observe  the  formation  (between  the 
contours  along  the  CM  within  the  values  at  t  =  20.2  hr,  Log.nT  =  G.  0  and  6.4)  of 
a  fast  reverse  MIfD  shock  wave.  See  als  >  D'Uston  et  al,  where  this  formation 
was  desci  ibed  in  more  detail.  At  t  =  40.2  hr,  the  shock  is  identified  between 
the  contours  of  Log^T  =5.6  and  6.0.  The  reverse  shock  reaches  1  AU  at  about 
t  =  60.  1  hr.  Note  that  the  forward  shock  becomes  attenuated  as  the  position  angle 
of  an  observer  increases  to  each  side  of  the  CM.  The  reverse  shock  is  much 
narrower;  it  too  becomes  attenuated  at  much  smaller  position  angles.  It  is  of 
great  importance  to  emphasize  the  substantial  amount  of  heliolongitudinal 


Bl.  D'Uston,  C.,  Dryer,  M. ,  Han,  S.  M. ,  and  Wu,  S.  T.  (1981)  Spatial  structure 
of  flare-associated  perturbations  in  the  solar  wind  simulated  by  a  two- 
dimensional  numerical  MHD  model,  J,  Geophys,  Res.  86:525-534. 
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Figure  B2.  Temperature  Contours,  Together  With  Superposition  of  Deformed 
IMF  (at  t  >  0  hr).  Following  Introduction  of  the  Shock  Wave,  Vs  -  3000  km 
soo'l,  At  -  3400  sec.  Numbers  are  logjQ  (°K) 


spreading  (beyond  90°)  accomplished  by  the  shock  whose  original  "angle  of  spread" 
at  18HS  was  only  24°.  It  is  also  of  great  interest  to  note  that  the  region  of  IMF 
polarity  reversal  (t,  -  40.2  hr)  takes  place  in  the  vicinity  of  the  reverse  shock. 

We  believe  that,  for  the  present  case,  the  reverse  shock  is  a  parallel  shock 
where  the  plasma  temperature  gradient  is  the  steepest,  hence  there  is  no  intro¬ 
duction  of  an  azimuthal  IMF  component.  At  tj>  <  30°  (that  is,  east  of  the  flare's 
C’M)  the  reverse  shock  is  quasi-parallel,  hence  an  azimuthal  IMF  is  required  on 
the  downstream  (anti-sunward)  side  of  the  shock.  The  subsequent  bending  of  the 
riUF  must  be  compatible  with  the  azimuthal  component  generated  by  the  eastward 
flank  of  the  fast  forwa rd  shock  (outside  the  computational  domain  at  t  ;  40.  2  hr). 
The  net  result  is  the  spatially-deformed  IMF  as  shown  near  (and  eastward  of)  the 
CM  in  Figure  B2(c). 

Before  leaving  the  temperature  contours,  we  wish  to  point  to  another  impor¬ 
tant  result:  the  outward  movement  near  CM  of  a  region  of  temperatures  that  are 
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cooler  than  the  pre-disturbed  values  at  the  same  points.  We  call  the  readers' 
attention  to  the  contours  that  intrude,  as  it  were,  sunward  from  the  local  mini¬ 
mum.  As  time  progresses,  these  sunward  contour  "intrusions"  become  more 
elongated.  It  can  be  seen,  by  comparison  with  the  contours  for  t  =  0  hr  [  Figure 
B2(a)l  that  the  temperatures  are  indeed  cooler  than  their  undisturbed  values 
along  this  narrow  heliolongitudinal  corridor. 

We  have  called  this  IMF  and  temperature  behavior  to  the  readers1  attention 

because  of  its  relevance  to  the  observations  of  magnetic  "clouds"  (Klein  and 
B2  B3 

Burlaga;  Burlaga  et  al.  ,  and  others).  We  will  return  to  this  point  in  several 
additional  discussions,  regarding  the  density,  total  pressure,  and  plasma  beta 
(fi  =  16 v  nkT/B2). 

Figure  B3  shows  the  contours  of  Log^n  (cm  2)  together  with  the  super-im¬ 
posed  IMF  from  Figure  Bl.  The  large  extent  of  the  low  density,  overexpanded 
regions  behind  the  fast  forward  MUD  shock  has  been  pointed  out  in  earlier  work 
(Wu,  Han,  and  Dryer;  4  Dryer,  Wu,  and  Han;  ^  and  D’l’ston  et  al.*4*)  and  is 
evident  here  again.  Here,  however,  we  can  relate  particularly-low  minima  with 
the  IMF  twists  and  turns  as  outlined  during  the  temperature  discussion. 

9 

We  also  call  attention  to  contours  of  plasma  beta  ((1  =  16*  nkt /B  )  in  Figur  e 
B4.  This  parameter  is  of  importance,  not  only  in  the  context  of  discussions  of 
the  magnetic  clouds  (Klein  and  Burlaga;  c  Burlaga  et  al.  ),  but  also  in  the  study 
of  the  structure  and  foreshock  char  acteristics  of  both  forward  and  rever  se  shocks. 
Plasma  f)  (ratio  of  thermal  to  magnetic  pressures)  is  also  of  great  importance  in 
the  eventual  delineation  of  a  hybrid  linkage  of  kinetic  and  MUD  modelling  of  inter  ¬ 
planetary  disturbances  in  a  self-consistent  way.  Figures  B4(b)  and  B4(c)  clearly 
show  the  attainment  of  extremelv-large  values  of  ,i  ->'*  1  that  are  co-spatial,  mure- 
or-less  with  the  IMF  reversa  1  regions.  We  also  set'  that  the  opposite  extreme, 

(1  <<  1,  occurs  within  the  much  larger  region  of  the  disturbed  flow. 


B2.  Klein,  I..W.,  and  Burlaga,  L.  F.  (1982)  Interplanetary  magnetic  clouds  at 
1  AF,  J.  Geophys.  ties.  87:613-624. 

B3.  Burlaga,  L.  F. ,  Sutler,  K. ,  Mariani,  F. ,  and  Schwenn,  K.  (1981)  Magnetic 
loop  behind  an  interplanetary  shock:  Voyager,  Helios,  and  IMP  8 
observations,  J.  Geophys.  Bes.  86:6673-6684. 

B4.  Wu,  S.  T. ,  Han,  S.  M.  ,  and  Dryer,  M.  (1979)  Two-dimensional,  time- 

dependent  MHD  description  of  interplanetary  disturbances:  simulation  of 
high  speed  solar  wind  interactions.  Planet,  Space  Sci.  27:255-264. 

B5.  Dryer,  M. ,  Wu,  S.  T. ,  and  Han,  S.M.  (1980)  Two-dimensional,  time- 

dependent  MHD  simulation  of  the  disturbed  solar  wind  due  to  representative 
flare-generated  and  coronal  hole-generated  disturbances,  Geofis.  Intern. 
19:1-15.  - 
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Figure  B3.  Density  Contours,  Together  With  Superposition  of  Deformed  IMF 
(at  t  >  0  hr),  Following  Introduction  of  the  Shock  Wave,  Vs  =  3000  km  sec-1, 
At  =  5400  sec.  Numbers  are  Log^Q  (particles  cm~3) 


As  discussed  in  several  earlier  papers  (Dryer,  Wu,  and  Han;®5  Wu,  Dryer, 
B6 

and  Han  ),  the  results  of  n,  T,  Ptotaj .  and  (3  appear  to  replicate  the  character¬ 
istics  of  magnetic  clouds  as  observed  and  delineated  by  Klein  and  Burlaga.®^ 

The  characteristics  to  be  discussed  are  shown  for  only  one  case:  V  -  3000  km 

-1  s 
sec  (At  =  5400  sec),  which  is  representative  of  the  entire  range  of  initial  shock 

velocities  mentioned  earlier  in  the  main  text.  It  should  be  remarked  that  we 

made  no  special  efforts  to  arrive  at  characteristics  discussed  by  Klein  and 

Burlaga,  thus  they  were  derived  in  a  completely  independent  way.  We,  therefore, 

will  continue  to  describe  our  results  with  comparative  remarks,  regarding  the 

magnetic  clouds,  left  for  a  discussion  later  in  this  appendix. 


B6.  Wu,  S.  T. ,  Dryer,  M. ,  and  Han,  S.  M.  (1982)  Non-planar  MHD  model  for 
solar  flare-generated  disturbances  in  the  heliospheric  equatorial  plane, 
submitted  to  Solar  Physics. 
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Figure  B5  shows  the  fractional  temperature  change,  (T  -  T  )/To>  as  the 
shock  and  post-shock  flow  pass  through  one  quadrant  of  the  computational  domain. 
With  reference  to  the  half-tone  scale  given  to  the  right  side  of  Figure  B5(a),  we 
note  that  a  very  high  Mach  number  shock  (which  is  not  unexpected  since  the  initial 
velocity  was  so  large)  is  seen  at  the  flare's  CM.  The  initial  heating,  of  course, 
gradually  decays  as  time  progresses;  eventually,  an  expanding  region  of  cooler 
gas  appears  as  seen  at  t  =  80  hr  in  Figure  B5(d).  This  result  can  also  be  seen 
in  the  contour  plot  of  Figure  B2(d)  which  shows,  at  t  =  80  hr,  temperatures  as 
low  as  10^ ‘  ®  K,  as  compared  w  ith  original  values  i  Figure  B2(a)l  of  ~  10^'  °K 
at  the  same  location  close  to  the  flare's  CM. 

Figure  B6  shows  the  fractional  density  change,  (n  -  ng)/n  corresponding 
to  the  same  times  shown  in  the  previous  figure.  Once  again,  as  with  Figures  B2 
and  B3,  we  can  compare  the  shock  compression  and  heating,  followed  by  adiabatic 
cooling  and  expansion  in  the  sunward  direction  until  the  reverse  shock  is 
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Figure  B5.  Fractional  Temperature  Change,  (T  -  T  )/T  ,  Within  Part  (0  s  <t  S 
00°)  of  the  Calculationai  Domain.  Note  that  the  central  p8rtion  of  the  fast  forward 
shock  has  left  the  domain.  .Also,  note  the  large  region  with  no  shading  (see  half¬ 
tone  scale  to  the  right  of  the  lower  right  panel)  with  significant  cooling  relative  to 
the  ambient,  undisturbed  temperature:  (T  -  T  )/T  -  -0.3 


encountered.  The  faster  velocity,  sunward  of  the  reverse  shock  (Figure  3),  flows 
through  this  discontinuity  where  a  much  lower  level  of  heating  and  compression 


(  a  I  ',Hf  ’U.’MBOO 


lb)  !;"t  ?  ./•■j."  r.-tv:'- 


‘»0K  .0810  1. ISO 


o>0M  .0810  j .  ;RL- 


( c )  *  :«f  60.  I*i0i, so  0EN51'1  1 0 '  00  - 1  i  Id )  "■">  8C.JH0U06  DEN811»  DOO-l- 


.0010  i.inu  FR0M  .0010  i.iou 

Figure  Rfi.  Fractional  Density  Change,  (n  -  n0)/n0.  Within  Part  (Oj  90°)  of 
the  Computational  Domain.  Note  that  a  substantial  constantly-expanding  region 
of  depressed  density  levels  follows  within  the  two  shocks'  wake.  Part  of  this 
region  coincides  with  the  depressed  temperature  region  I  compare  Figure  R(i(d) 
with  Figure  R.">(d)l 


takes  place.  The  overexpansion  that  follows  behind  the  reverse  shock  (and  its 
attenuated  longitudinal  extent)  is  apparent  in  each  of  the  boxes  in  f  igure  Rfi  (as 
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noted  by  the  negative  values  in  the  keyed  half-tone  scale  next  to  each  box).  Note 
that  the  regions  of  depressed  (lower  than  original  steady-state  values)  tempera¬ 
tures  and  densities  [  Figures  B5(d)  and  RG(d)l  overlap  in  a  scale  size  that,  at  1 
AU.  is  about  0.25  AU  in  longitudinal  extent.  Comparison  of  these  two  figures 
shows  that  the  length-to-width  ratio  of  the  depressed  density  region  is  oriented 
in  the  ^-direction,  whereas  that  of  the  depressed  temperature  region  is  oriented 
in  the  r -direction.  We  suspect  that  the  region  of  overlapping,  together  with  part 
of  the  region  of  OUF  reversals,  may  constitute  a  magnetic  cloud  whose  plane  of 

minimum  variance  (of  the  IMF  component  perpendicular  to  it)  is  reasonably  close 

B2 

to  our  equatorial  plane.  Klein  and  Burlaga  noted  that  their  set  of  70  candidate 
magnetic  cloud  events  included  some  (~  25)  for  which  the  minimum  variance 
direction  was,  indeed,  nearly  normal  to  the  ecliptic  plane. 

Figure  B7  shows  the  fractional  change  of  total  pressure,  (P  ,  ,  ~  P.  ,  ,  )/ 

2  tot  3.1  totaiQ 

P(otai  .  where  Ptotaj  2nkT  +  B  /8 There  is,  of  course,  a  rapid,  precipitous 
increase  as  a  result  of  the  fast  forward  shock  compression  and  heating.  An 
observer  at  1  AU  along  or  near  the  CM  would  observe  this  rapid  rise,  followed 
by  a  gradual  decrease  over  a  period  of  about  3  days,  as  suggested  by  the  results 
in  Figures  B7(a)  -  B7(c),  Thereafter,  for-  about  2  days,  the  total  pressure  would 
be  predicted  to  decrease  below  its  steady-state  value.  This  result  appears  to 
contradict  earlier  discussions  of  cloud  expansion  due  to  high  total  pressure  within 
them.  In  fact,  it  does.  The  earlier  discussions  implicitly  assume  that  static 
equilibrium  is  pr  esent.  By  virtue  of  the  time  variance  of  all  parameters  at  an 
observing  spacecraft  (see  Klein  and  Burlaga  c),  the  flow  is  time  variant  in  the 
spacecraft  frame.  Hence,  a  simplistic  determination  of  a  gradient  of  total  pres¬ 
sure  is  insufficient  to  arrive  at  a  valid  conclusion  concerning  "cloud"  expansion, 
cooling,  etc.  Here,  we  have  assumed,  implicitly,  that  reconnection  need  not 
take  place;  if  true,  extraneous  bubbles  might,  in  principle,  float  with  the  solar 
wind  or  even  move  relative  to  the  plasma  like  a  projectile.  W c  think  it  premature 
at  this  time  to  require  such  exotic  systems.  Thus,  we  believe  that  a  simple  time- 
variant  system  as  described  here  is  sufficient  for  a  fundamental,  self-consistent 

understanding  of  the  basic  physical  process  for,  at  least,  the  shock-associated 

B  2 

magnetic  clouds  described  by  Klein  and  Burlaga. 

The  discussion  of  temporal  phenomena,  as  found  by  the  MHD  model,  would 

be  incomplete  without  a  presentation  of  the  plasma  and  field  parameters  at  vari- 

B  6 

ous  observational  points.  In  the  paper  by  Wu,  Dryer,  and  Han,  we  used  the 
example  of  Vg  =  1000  km  sec  1  to  present  all  parameters  (including  velocity  and 
magnetic  fields)  at  various  points,  each  one  considered  separately.  Here,  we 
use  a  variation  of  this  former-  presentation  in  order  to  emphasize  heliolongitudi¬ 
nal  attenuation  of  the  basic  physical  quantities:  radial  velocity,  density,  tempera¬ 
ture,  and  total  magnetic  field  at  0.3  as  well  as  at  1.0  AU. 
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Figure  B8.  Time  Series  of  Radial  Solar  Wind  Velocity,  Density,  Temper¬ 
ature,  and  Total  IMF  Magnitude  at  0.3  All.  (a)  -  (d)  CM  of  flare-generated 
shock,  (e)  -  (h)  33°W  of  CM,  and  (i)  -  (1)  90°W  of  CM  (Contd) 
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(I'rurn  a  practical  point  -of-viow,  the  two  Helios  spacecraft  move  very  quickly 
through  their  0.3  AV  perihelia,  lienee  their  hypothetical  counterpart  would 
"sample",  only  a  short  portion  of  the  200-hr  series  shown  in  Figure  B7.  ) 

Figures  B8(n)  -  B8(d)  show  the  radial  velocity,  density,  temperature,  and 
total  magnetic  field  at  the  CM  of  the  canonical  flare  that  produced  a  3000 
km  sec’'  shock  wave  at  18Rs ,  as  used  in  the  present  simulation.  The  out¬ 
standing  featur  es  are:  (1)  prompt  arrival  about  4  hr  later  of  the  fast  forward 
shock  followed  at  t  ~  10  hr-  by  the  initiation  and  build-up  of  the  fast  reverse  shock; 
(2)  rapid  compression,  followed  by  an  extensive  rarefaction  to  greatly-reduced 

densities  for  about  50  hr  before  the  recompression  and  very  gradual  return  to  the 

6 

prr-disturbed  values  of  density;  (3)  rapid  heating  to  a  temperature  of  ~  24  <  10 
°K,  which  is  common  in  a  large  flare;  and  (4)  magnetic  compression  and  rare¬ 
faction,  accompanied  by  the  strongly -deformed  topology  shown  in  Figure  Bl,  with 
a  recovery  to  only  one-half  its  pre-disturbance  value  at  the  termination  of  the 
simulation.  This  is  a  real  result  and  not  a  numerical  artifact;  the  azimuthal  and 
meridional  IMF  components  return  very  nearly  to  their  pre-disturbance  values 
after  the  strong  shock  perturbation,  but  the  radial  component  appears  to  respond 
in  a  wiggly  fashion  as  the  slow  mode  wave  moves  outward  as  well  as  in  the  longi¬ 
tudinal  direction  as  discussed  earlier  in  the  discussion  of  Figures  3(e)  -  3(g). 

Panels  (e  -  h)  and  panels  (i  -  1)  of  Figure  B8  show  the  same  parame¬ 
ters  at  3 3 °\V  and  90°W,  respectively,  of  CM.  Noting  the  increase  of  the  vertical 
scale  [  especially  for  the  case  of  90°\V  in  Figures  B8(i)  -  338(1)1 ,  the  delay  (of 
arrival  time)  and  attenuation  of  the  peak  velocity,  density,  and  temperature  at 
the  forward  shock  (even  at  90°W)  arc  clearly  seen.  Note,  however,  that  the  peak 
IMF  is  increased  at  33°W  (to  G10  nT)  relative  to  its  peak  value  at  CM  (385  nT). 
The  reason  for  this  result  is  straightforward  when  one  keeps  in  mind  the  fact  that 
the  shock  propagates  asymmetrically  in  a  simultaneous  radial  and  heliolongitudi¬ 
nal  fashion.  At  CM,  then,  the  shock  is  very  nearly  quasi-parallel  -  but  not  quite 
so;  thus  the  IMF  jump  is  [  from  Figure  B8(d)l  only  from  235  nT  to  385  nT.  At 
33  W,  however,  the  shock  normal  inclination  to  the  IMF  (  0n  g )  is  larger  than  at 
CM;  hence  the  magnitude  jump  I  Figure  B8(h)l  is  larger,  Trom  235  nT  to  810  nT 
despite  some  attenuation  of  the  shock  velocity  in  that  direction.  By  90°W,  how¬ 
ever,  the  shock  velocity's  attenuation  overcomes  the  influence  of  a  quasi-perpen¬ 
dicular  value  of  g,  and  the  jump  [  Figure  B8(l)l  is  from  235  nT  to  only  297  nT. 

Further,  wc  explored  a  simple  variation  to  our  ansatz  problem.  What  would 
be  the  Karth-oriented  effect  on  c  and  VBg  if  the  Bg  component  at  18Kg  were 
initially  an  order  of  magnitude  larger  (10  nT)  than  the  value,  Bg  1  nT,  as  used 
heretofore  (Table  1)?  To  answer  this  question,  wo  performed  another  calculation 
(about  3  min  on  the  Cray-1,  33  min  on  the  Cyber-750)  with  this  change  to  the 
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parameter -s  in  1  able  1.  The  results  an-  shown  in  F igure  B9.  Contrary  to  a 
simolistic  kinematic  expectation,  we  did  not  obtain  a  larger  e  simply  as  a  result 
of  an  inilially-larger  southward  IMF  component.  Figure  B9(a)  shows  that  the 

19  _  1 

magnitudes  of  the  two  maxima  arc  reduced  below  the  10  '  erg  sec  threshold 
that  had  otherwise  been  breached  in  Figure  ti(a).  On  the  other  hand,  VB^  was 
increased  by  more  than  an  order  of  magnitude,  which  would  be  in  agreement  with 
a  kinematic  interpretation  that  is  not  dependent  upon  a  realistic  evolution  and 
kinking  of  the  IMF. 


Figure  110.  ener  gy  Flux  Parameter,  e,  and  the  Product  VB^  at  1  AU  For  the 
Case:  Vs  3000  km  sec"l,  At  -  3400  sec,  Bgg  =  10  nT  (Instead  of  1  nT,  as 
Fsed  Karlier).  (a)  CM,  (b)  33°W  of  CM,  and  (c)  90°W  of  CM 

This  latter  point  concerning  the  evolution  of  the  IMF  brings  us  to  our  final 
comments  regar  ding  the  earlier  "magnetic  cloudl  discussion.  We  cross-plotted 
the  IMF  components  against  each  other  (without  performing  a  minimum  variance 
computation  in  this  early  exploratory  work)  for  the  case  Vg  =  3000  km  sec"*, 
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10  ti  l'.  The  results  ;tri-  shown  in  I  inure  H10  in  the  form  of 


At  .400  sir;  li() 

polar  plots  usstiimnc  that  tin-  equatorial  plane  is  a  reasonable  approximation  to 
tile  true  plane  of  minimum  variation.  It  earl  be  seen  that,  indeed,  the  excursion 
of  Ii^j  1 1  o e s  not  exceed  ±  20  nl'.  whereas  the  11^  component  has  an  excursion  to 
nearly  80  nl'  in  the  eastward  direction.  It  is  of  interest  to  note  that  a  rotation  of 
nearly  !<0°  of  the  IMF  field,  projected  into  the  equatorial  plane,  takes  place.  (It 
is  also  seen  in  l-  igure  151(c),  at  positions  east  of  CM  for  the  previous  case,  that 
a  180°  rotation  takes  place.) 


I'M)  SIMULATION  AT  1.0  AU 
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Figure  B10.  Polar  Plots  of  IMF  at  1  AC.  CM,  For  the  Case: 

3000  km  sec'*,  At  -  3400  sec,  B#  -  10  nT.  Note  the  scale  change 
between  B,  and  B() 
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